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Abstract 

We report a novel practical method to determine exchange parameters 
by using experimental susceptibility data in a relatively narrow temper- 
ature region and a low order high-temperature-expansion equation. This 
method is applied to a square lattice and a CaV40g-type lattice, and its 
accuracy is discussed. 
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When we describe magnetic properties of a matter by a Heisenberg model, we 
must determine exchange parameters in the model to fit the matter. It is some- 
times performed by the high temperature expansion (HTE) of the susceptibility 
of the Heisenberg model and by the comparison of it to the experimental data. 
However meaningful experimental data are often restricted in a low temperature 
region because of, e.g., the temperature dependence of the lattice constant or 
structure. Further it becomes hard to calculate high order terms of the HTE 
because of striking increase of the number of diagrams if the model includes 
many exchange parameters. In some of such cases, an ordinary method of the 
direct fitting may not estimate exchange parameters with sufficient accuracy. 

In this article, we report a novel practical method to determine exchange 
parameters by using experimental susceptibility data in a relatively narrow tem- 
perature region and a low order HTE equation. The essence of this method is 
that we construct an experimental formula for the experimental susceptibility 
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data in a power series of 1/T and then compare the coefficients of this formula 
to those of the HTE equation derived from the assumed Heisenberg model. 

We denote the experimental data for the magnetic susceptibility by x EXP {T) 
and the experimental formula in a power series of 1/T by x{T). To make ex- 
pansion coefficients dimensionless we use x = Tq/T as the expansion parameter, 
where T is a constant temperature chosen arbitrarily. Then x(T) is written as 

X (T) = lim x (n) {T), 

p n 

X {n \T) = -Y,A m x m , (1) 

where C is Curie constant determined by experiments. The zeroth expansion 
coefficient is set as Ao = 1 and the others (A m 's for m > 1) are fitting parameters 
to be determined. 

Using the same coefficients A m 's, we define the following quantity: 

oo 

<f> m (x) = A m +J2A m+ ix l (2) 
l=i 

for each m. In the high temperature limit, it reduces to a coefficient as 

<f> m {0) = A m . (3) 
This quantity <fi m (x) is also defined by the recursion equation: 

T 

M x ) = q x(t), 

<t> m {x) = ((/> m -i(x) - A m -i)- (4) 

X 

for positive integer m. 

The experimental data corresponding to 4> m (x) is similarly defined in the 
same recursive transformation: 

tf xp {*) - ^x EXP (n 

4>l xp (x) = (^f(x)-A m _ 1 )l. (5) 

x 

The quantity 4>m XP 1S actually determined only when A m _i is known. It is also 
noticed that 4> E ^ P is defined only at experimental data points in {x}; e. g. it 
is not defined at very high temperatures (a; ~ 0). <f> m must fit <^ , since x 
is constructed to fit x EXP an d 4>m is produced by the same transformation as 

<b EXP 

In the following, we inductively determine {A m } together with {(j) 1 ^ p (x)} . 
For m = 0, 4>q XP (x) is defined by eq. ([)]) and A — 1 is a definition. Assuming 
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that <fi m -i(x) and A m -i are known, we construct (f> 1 ^ l xp (x) and A m . (f>^ l xp (x) 
is simply given by eq. (0). To determine A m , we consider the high temperature 
limit (x — » 0) for 4>f n (x) which corresponds to eq. (||) for 4> rn (x). Since 
4 l m XP ( x ) does not have data points in the vicinity of x=0, we extrapolate the 
graph of (f)f J ^ lp {x) to x = by using a polynomial fitting function: 

L 

fm(%) = y^QWz': (6) 
1=0 

where coefficients a m /'s (Z = 0, • • • , L) are fitting parameters to be determined 
for each m. The order L of the polynomial is empirically chosen as 8 to 10 and 
is inspected with some examples as mentioned later. A m is then determined 
by A m ~ / m (0) = a m o corresponding to eq. Thus we have completed the 
inductive definition and have constructed an experimental HTE formula (lb). 

On the other hand, we start from a Heisenberg Hamiltonian H which is 
assumed to describe the magnetic properties of the material. The assumed 
Hamiltonian includes a set of exchange parameters { Ji}. For this Hamiltonian, 
we calculate the theoretical HTE formula 

X HTE (T) = % £ F m x™ (7) 
m=0 

by the standard diagrammatic method, jlj The expansion coefficients -F m 's are 
functions of the set of exchange parameters {Ji}- 

Lastly we compare the theoretical HTE formula (Q) to the experimental 
formula (|l|) . If the assumed Hamiltonian exactly describes the material, both the 
coefficients must be the same. Hence we can determine the exchange parameters 
by solving the set of equations {F m {{Ji}) = A m }. In real cases, however, the 
model Hamiltonian may not perfectly describe the matter and the experimental 
data may include some errors. We usually must be satisfied with a set {Ji} 
which gives approximate matching between {F m } and {A TO }. 

We inspect the method mentioned above for the Heisenberg model in a simple 
square lattice. Only an exchange parameter J is between nearest neighbor spins. 
As data for susceptibility, we use the result of Quantum Monte Carlo (QMC) 
simulation by Troyer et al. M instead of data of a real experiment. We choose 
the expansion parameter as x = 1/T in the unit of To = 1 and use 12 QMC 
data points in 0.05 <T< 1.5. The data is for J = 1 and C = 0.25. 

Following the above procedure, we obtained coefficients of HTE as Ai = 
-0.994, A 2 = 0.483 and A 3 = -0.135 for L = 8. These values are little changed 
when the order L of the fitting function changes from 6 to 10. We show <p^ l xp (x) 
along with f rn (x) for m = 1 to 3 in Fig. [I]. 

The corresponding theoretical coefficients are calculated as F\ = — J, F 2 = 
J 2 /2 and F3 = — J 3 /6. In this one-parameter case, equation F\(J) = A\ is 
sufficient to determine J. The solution is J = 0.994, which is close to the 
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Figure 1: Transformed QMC data (f>^ l xp (x) (circle) and fitting function f m (x) 
(line) with m = 1 to 3 for a square lattice. 

exact value J — 1. We estimate the accuracy of A m by the deviation d m = 
2\(A m - F m (l))/(A m + F m (l))\. It is given as d rn = 0.006, 0.017 and 0.070 for 
m = 1, 2 and 3, respectively. This result, i. e. this method, is fairly accurate 
despite the small number of the data points. 

We also carried out the same analysis by restricting the temperature region 
as 0.05 < T < 0.9, which includes only 9 data points. This region is below 
the maximum point of x (T) and is rather a low temperature region. The 
present method gives A x = -1.019, A 2 = 0.541 and A 3 = -0.185 for L = 8. 
The corresponding deviations are d m — 0.019, 0.079 and 0.327 for m = 1, 2 and 
3, respectively. This estimation is practically sufficient to determine J from A\. 

Next we apply our method to a Heisenberg model on a lattice which is a 
simplified two-dimensional model used for CaA^Og. This model includes two 
exchange parameters, Jo for the plaquette link and J\ for the dimer link. Theo- 
retical coefficients are obtained as J\(Jo, Ji) = — (2J + Ji)/4 and F 2 (J , Ji) = 
(4Jo — Ji)Ji/16. Instead of experimental data, we use QMC data @ and ap- 
ply our method to determine Jo and J\. We used 15 QMC data points in 
0.075 < T < 1.5 and the data are for C = 0.25 in the unit of T = 1. We 
examine two typical cases of Jo/Ji =1.0 and Jo/Ji = 0.5. The system has no 
spin gap in the former case, while it has in the latter. Q 

For the data of (Jo, Ji) = (1, 1), our method produces coefficients, A\ = 
—0.750 and A 2 = 0.166. The fitting was done by the 8th order polynomial 
(L = 8), / m 's. The corresponding exact values are F\ = —0.750 and F 2 = 0.188, 
and the deviations are d\ = 0.00 and d 2 — 0.12. The complete matching between 
A\ and F\ over 3 digits are occasional, since the last digit changes if we change 
L from 8 to 9 or 10. We calculated the exchange parameters as (Jo, Ji) = (0.83, 
1.34) or (1.17, 0.66) by solving J\(J , Ji) = A x and F 2 (J , J t ) = A 2 . 

Similarly we applied our method for the data of (Jo, Ji) = (0.5, 1). We 
obtained HTE coefficients as A\ = —0.500 and A 2 = 0.052. The corresponding 
exact values are F\ = —0.500 and F 2 — 0.063, and the deviations are d\ = 0.00 
and d 2 ~ 0.19. The obtained exchange parameters are (Jo, Ji) = (0.46, 1.08). 

The estimation gives better result for a square lattice than for a CaA^Og 
lattice. A reason may be small HTE coefficients included in the latter lattice. 
Such small coefficients seems to appear, when the coordination number is small 
so that the spin fluctuation is strong. Large fluctuation is also a cause of the 
spin-gap formation. The present method is applicable to gapful systems as well 
as gapless systems as long as we only need coefficients of low orders. 

We thank K. Kubo for useful discussion on high temperature expansion 
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4 



This work is partially supported by the Grant-in-Aid for Scientific Research 
from the Ministry of Education, Science and Culture, Japan. The computation 
in this work was carried out partially by using facilities of the Supercomputer 
Center, Institute for Solid State Physics, University of Tokyo. 

References 

[1] For example, G. S. Rushbrooke, G. A. Baker, Jr. and P. J. Wood: Phase 
Transition and Critical Phenomena, eds. C. Domb and M. S. Green (Aca- 
demic Press, London and New York, 1974) Vol. 3, Chap. 5, p245. 

[2] M. Troycr, H. Kontani and K. Ucda: Phys. Rev. Lett. 76 (1996) 3822. 



5 



